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A discretisation scheme that preserves topological features of a physical problem is extended so 
that differential geometric structures can be approximated in a consistent way thus giving access to 
the study of physical systems which are not solely topological theories. Issues of convergence and 
a numerical implementation are discussed. The follow-up article covers the resulting discretisation 
scheme with metric data. 
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Problems in physics and engineering can often be formulated using the coordinate independent language of differ- 
' ential forms I*!. The setup can be characterized in terms of the following objects: a manifold M on which p-forms 
(anti-symmetric tensor fields) ujp are defined, the presence of a differential operator d that maps a p-form to a (p+ 1)- 
form, a wedge product. A, which allows a p-form and g-form to be multiplied to produce a (p + (7)-form in a way that 
04 allows a Leibniz type rule to be valid for the operator d when it acts on the product and the Hodge star operator 
(-H . The Hodge star operator requires that the manifold M have a metric. It maps a p-form to a {D — p)-form where 
^ D is the dimension of the manifold M. Using the ★ operator a scalar product between p-forms can be introduced. 
This immediately allows the adjoint d'^ of d to be defined. The geometric discretisation scheme 0, 0, S IE IE 13 

■ constructs an analogue model involving discrete objects for the continuum system described. Each of the elements 
(wp, d, A, d^) now have discrete analogues. 

Two maps allow the approximation scheme to be defined precisely. The de Rham map that maps continuum vari- 
, ^ \ ables, such as Wp, to discrete variables Cp "chains" and the Whitney map which maps discrete objects cjp to continuous 
. variables W^[a'p] "cochains". Results establishing the nature of the approximation and questions of convergence have 
' been established 5] . One generic feature of the scheme needs to be stressed. The discrete analogue of the Hodge star 
, mapping CTp *(Tp produces a different space of objects. Hence in order to properly work out discrete analogies of 

■ d'l' it becomes necessary to work in the space which is symbolically Mk ® *Mk where Mk is the discretisation of the 
, manifold M. There is a doubling of spaces introduced 

' In previous work it has been shown that the geometric discretisation scheme captured topological features of a 
\ system very nicely l9l. Extension of the ideas to accommodate fermions, in a Dirac-Kahler framework, have also 
I been carried out |lOl Hil l . Here we turn to the practical problem of geometrical information present in the original 
O i' continuum problem in the discrete approximation scheme described without spoiling topological structures that it 
O , already encodes. The approach we proceed to introduce, tackles the problem of defining what is known as a discrete 
' exterior calculus and different attempts to tackle parts of this problem are available, Bossavit has introduced a 
^ discrete analogue of contraction by considering the dual mapping "extrusion" in the context of the electromagnetic 
theory. Hirani and co-workers [l^ have developed part of a discrete exterior calculus which is closely linked in spirit 
to the topological theory (the geometric discretisation described in I. being an example), but they do not study the 
' Lie derivative which is of much importance here. 
■ ■ ' In contrast to other works, in our approach the Whitney map which maps discrete objects to continuum variables 
plays an essential role. A definition of a product space of chains allow the notion of contraction to be defined which 
is essential for introducing analogues of the Lie derivative, an operation that describes the way objects change. The 
structure constructed gives rise to a consistent discrete exterior calculus, which is to our knowledge a new result. 



I. THE TOPOLOGICAL DISCRETISATION SCHEME 



The original geometric discretisation scheme captured topological data. For the simplical complex, it was developed 
by several authors 0jlESIES0@- Here we will consider the hyper-cubic version of it 0| as our starting point. 
An alternative derivation of the Whitney elements is introduced which suits our purpose. First, we point out the 
important shortcomings of the original geometric discretisation method: 

• The discrete wedge was non-associative. 
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• Operators are only exact as operations on the discrete chains, (i.e there is no direct correspondence or "functor" 
between the discrete operations and the continuum ones). 

• Whitney elements play a spectator role and so the theory is topological only. 

Recently, Samik Sen [l^ has shown that by replacing a simplicial complex by a hyper-cubic complex, two new features 
emerge compared to the original proposal: the Hodge star operator satisfy ★^★^ = Id and second, the associated 
dual complex being hyper-cubic as well, it can be provided with explicit expressions for the Whitney elements. 

The second item of the list above is crucial to us, because we are going to take the Whitney forms as the basic 
building blocks of the theory. Then, we use the product space of chains to generate the appropriate collection of forms 
which will in turn allow us to consider a discretised system with symplectic geometric data. 

In contrast, the original scheme does not utilize the Whitney forms beyond providing a space that gives rise to the de 
Rham co-homology complex, so we need not know the explicit form of the Whitney elements to proceed topologically. 

Here, we put the Whitney elements and the various combinatorial manipulations of them at the forefront of the 
formulation. Here, a discrete operation is exact if when mapped to the continuum its algebraic properties hold exactly. 
Then in turn, by accommodating the other operations in the theory (A, d, ★) we can show that the combinatorics of 
Whitney elements that we introduce is such that the topological data is still captured. We find that: 

• The basis of Whitney forms encodes the topology and is very constraining that way, but allows for the product 
space of chains. 

• Secondly, the problem of the non- associativity of the wedge is by-passed by not relying on the wedge but on a 
new product in order to define the contraction. 

These arguments both provide a strong motivation for considering the product space of chains. Before we turn to 
that, we review the geometric discretisation scheme. 

A. Differential forms, manifolds and complexes 

The general program of the topological discretisation is described at length in 9], On a manifold, we have fields 
(j}(p) which are differential forms (antisymmetric tensor fields), and the various operations of interest are: to consider 
a theory defined on an arbitrary manifold M of dimension D. Suppose the theory has fields (j>^{x), where x ^ M and 
p = 0, 1, . . . , D. The theory is constructed using the following objects which are defined on the manifold M : p- forms 
which are generalized antisymmetric tensor fields, the exterior derivative d : 4)^ (j)^^^, the Hodge star operator 
* : (pP ^ , which is required to define scalar products, and the wedge operator (pP A (f)'^ = (p^^^ . On a manifold 

M of dimension D, the operation (A, *, d) on p-forms, 4>^{p = 0, . . . , D), satisfy the following: 

1. (/.f A = {-If^" A 

2. d{4f A 09) = d4f A + {~lY(t>P A d4>'^. 

3. *(t)P = (fy^-P. 

4. = (-l)^P+i. 

5. d^ = 0, (d*)2 = 0. 

6. = {-l)D{p+mi ^d*. (dt is the adjoint of d). 
The following definitions are also important: 

• The Laplacian on p-forms Ap = dp-id* + dp_^^idp. 

• The inner product {(j)P, cfyP') = (f)P A *(j)P' . 

We would like to discretisc the fields of 4)P , the inner product (•, •) and the operators A, *, d such that discrete analogies 
of the above interrelationships hold. To do this it is necessary to first introduce a few standard ideas. We first discretise 
the manifold M by replacing M by a collection of discrete objects, known as simplices, glued together For p > 0, 
a p-simplex crtP) is defined to be the convex hull in some Euclidean space of a set of p + 1 points wq, wi , . . . , £ 
where the vertices are independent in that 
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FIG. 1: Dual complex L. 



implies that = for i — 0, . . . ,p where Aq, Ai, . . . , Ap are real numbers. Geometrically the barycenter of a given 
n-simplex cr*^"^ = {vo,vi, . . . ,Vn) is defined by 

n + I 

Thus the barycenter of cr^^^ = [vq, f i] is the midpoint between the vertices vq and vi. We can now describe a particular 

way that a given manifold M can be discretised. Let S" be a collection of simpiices {cj"''}, n = 0,1, . . . , D, with the 
property that the faces of the simpiices which belong to S also belong to it. The elements of S glued together in the 
following way is known as a simplicial complex|l7j , |18| | : 

1. CT,-"'' n ct'-'"'' = if uf"'', cr^''^ have no common face. 

2. cTj-"'' n aj'^^ 7^ if cTj-"'', crj'^'' have precisely one face in common, along which they are glued together. 

In many cases of interest (including all differentiable manifolds ^^), M can be replaced by a complex K which it 
is topologically equivalent to. K is then said to be a triangulation of M(Note this triangulation is not unique). In 
this way of discretising M, the building blocks are zero, one, . . . , D-dimensional objects, all of which are simpiices 
e.g. generalized oriented tetrahedra. 

The same manifold can be discretised using hyper-cubes or a CW complex. We will adopt the hyper-cubic approach 
which is the neater one, for the following reason which has to do with the construction of the complex dual L. First 
note that in contrast with FIG.l, the analogue construction for the dual complex in the case of the hyper-cubic 
complex leads to two superimposed complexes shifted by a half edges length, and in turn the Whitney elements are 
readily available [T^. 



B. The discrete operations. 



Let us first recall the construction of the discrete Hodge star for the simplicial complex which leads to L. In the way 
analogue to the continuum property (3. above), a map from a {D — p) dimensional object to a p dimensional object 
was introduced. Now consider the barycentres of the building blocks of the system K. We first construct {D — p) 
dimensional objects whose vertices are barycentres of a sequence of successively higher dimensional simpiices, where 
each simplex is a face of the following one. In other words {D —p) dimensional objects of the form {ap, CTp+i, . . . , ctd}, 
where (j„ is a face of Un+i- The orientation of these are set so as to be compatible with the manifold. Joining these 
objects together gives us the dual of a p. 

We can now give the general rule for mapping an n-simplex (t„ = [dq, ■ ■ • ^Vn] to a (Z? — n) dimensional object ( 
(D-n) cell) as follows: We think of f7„ as an element of a simplicial complex K. We have 

*K ■ bo, ■ • ■ , Vn] U[(f„, (T„+l, . • . , ctd], 

where (T„+i is the barycenter of an l)-simplex which has cr„ as a face, an+2 is the barycenter of an (n-|-2)-simplex 
which has an+i as a face and so on. These objects have to be coherently oriented with respect to [vq, ■ ■ ■ ,Vn\- The 
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set of these cells constitutes the dual space K of K. By this procedure, a discrete version of the Hodge star operation 
was constructed. Let us explain. The Hodge star operator involves forms. It maps p-forms in D dimensions to a 
[D — p)-form. The *k map involves not forms but geometrical objects. There is a simple correspondence relation 
between these two cases. Given a p-form, 0^, and a p-dimensional geometrical space, Ep, the p-form can be integrated 
over Ep to give a number. Thus Ep and (j)p are objects that can be paired. We can write this as a pairing 




In order to proceed, we need to introduce some more structure. We start by associating with a simplicial complex K, 
containing {o-p} {i = 1, . . . , Kp;p = 0, . . . , D) a vector space consisting of finite linear combinations over the reals of the 
p-simplices it contains. This vector space is known as the space of p-chains,Cp(-?ir). For two elements (7*, cr^ G Cp{K), 
a scalar product (ap,<j^) = (5] can be introduced. An oriented p-simplex changes sign under a change of orientation 
i.e. if (Tp = [vo, . . . ,Vp] and r is a permutation of the indices [0, . . . then [iV(o) , ■ ■ ■ , ^t(p)] = (~l)^bo) • • • ) Vp], with 
T denoting the number of transpositions needed to bring [^^-(o) , • • • , ^'r(p)] to the order [vq, . . . ,Vp]. 
Given the vector space Cp{K), the boundary operator can be defined as 

: Cp{K) ^ Cp-i{K). 

It is the linear operator which maps an oriented p simplex cr^P^ to the sum of its (p-1) faces with orientation induced 
by the orientation of tr^. li = [vq, . . . , Vp], then 

p 

daP = J2i-mvo,...,Vi,...Vp], 

i=0 

where [vo, ■ ■ ■ ,Vi, . . . ,Vp] means that the vertex Vi has been omitted from to produce the face "opposite" to it. 

Given that Cp{K) is a vector space, it is possible to define a dual vector space C'p{K), consisting of dual objects 
known as cochains; that is we can take an element of Cp{K) and an element C'p{K) to form a real number. Since the 
space Cp{K) has a scalar product, namely if (Jp.a^p e Cp{K) then {ap,a^) = Sij. We can use the scalar product to 
identify Cp{K) = Cp{K), so that we can consider oriented p-simplices as elements of Cp{K) as well as Cp{K). We 
can write our boundary operation as 

{[vo, . . . ,v^, . . . ,Vp],d^[vQ, . . . ,Vp\) = (-1)'. 
This suggests introducing the adjoint operation cLk defined as 

(rf^[uO, ...,Vi,..., Vp], [vo, ...,Vp]) = {[vo, ...,Vi,..., Vp], d^[vo, • • • , Vp]). 

This is the coboundary operator which maps Cp{K) — > Cp+i{K). 
Indeed we have 

d^[vo,...,Vp\ = '^[v,vo,...,Vp], 

V 

where the sum is over all vertices v such that [v,Vo, ■ ■ ■ ,Vp] is a (p+1) simplex. 

The boundary operators dx and the coboundary operator dx have the property DkQk = dxdx = 0. Furthermore, 

dn '■ Cp — » Cp+i, 
Ok ■ Cp ^ Cp-i. 

These operators are the discrete analogies of the operators d and (— * = which act on forms. 

These operators could be defined only when a scalar product ("metric") was introduced in the vector space Cp's. 
At this stage we have a discrete geometrical analogue of d, d'^ and We have also commented on the fact that the 
operation ★ maps simplices into dual cells i.e. not simplices. If the original simplicial system is described in terms of 
the union of the vector spaces of all p-chains then the space into which elements of the vector space are mapped by 
★ is not contained within this space, unlike the situation for the Hodge star operation on forms. It is found that this 
leads inevitably to a doubling of the fields when discretisation, preserving topological structures, is attempted. 
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C. Relation to the Whitney map 

We now need a way to relate a p-chain to a p-form. This together with a construction which hnearly maps p-forms 
to p-simphces allows us to translate expressions in the continuum to a corresponding discrete geometrical objects. We 
start with the construction of the linear maps from p-chains to p-forms due to Whitney 0] . In order to define this 
map, we need to introduce barycentric coordinates associated with a given p-simplex cjP . Regarding a'^ as an element 
of some R^, we introduce a set of real numbers (/io, • ■ • , IJ-p) with the property 

> 0, 

i 

A point X can be written in terms of the vertices of and these real numbers as 

p 

Note if any set oi Hi = then the vector x lies on a face of a^. One can think of x as the position of the center of 
mass of a collection of masses (/uo, . . . , fip) located on the vertices {vq. . . . ,Vp) respectively. Setting /ii = for instance 
means the center of mass will be in the face opposite the vertex Vi. The Whitney map can now be defined. We have 

: CP (if) ^ $P(AO, 

where ^^{K) is a p-form. If e Cp{K) then 

p 

W[aP] =p!^(-l)V,d^o A...dtiiA. ..dfip, 

i=0 

where dfj,i means this term is missing, and (/io, . . . , /ip) are the barycentric coordinate functions of a^. 

We next construct the linear map from p-forms to p-chains. This is known as the dc Rham map. We have 

: ^P{K) CP{K), 

defined by 

for each oriented p-simplex e K. 

A discrete version of the wedge product can also be defined using the Whitney and de Rham maps such that : 
CP{K) X Ci{K) CP+'i{K) as follows: 

xA^y^ A^{W^{x) A W^{y)). 

It has many of the properties of the continuous wedge product in that it is skew symmetric and obeys the Leibniz 
rule but it is non-associative. 

At this stage we have introduced all the building blocks necessary to discretise a system preserving topological 
structures. We summar ize t he properties of the maps introduced in the form of a theorem. We have 

Theorem (Whitney (l^ ) 

1. A^W'^ ^ Identity. 

2. dW^ = W^d^, where d : (t)P ^ <jf+\ 

3. L| W^^(a) =< a,l3>,a,l3(£ K 



4. d^A^ = A'^d. 
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This theorem shows how d can be considered as the discrete analogue of d. We now show how * can be considered 
as discrete analogue of *. For this we need barycentric subdivision. 

In order to construct the star map, two geometrically distinct spaces were introduced. The original simplicial 
decomposition K with its associated set of p-chains Cp{K) and the dual cell decomposition L with its associated set 
of p-chains Cp{L). These spaces are distinct. However both belong to the first barycentric subdivision of K. This 
allows the use the *^ operation if wc think of K and L as elements of BK. 

We proceed as follows. Let BK and L denote the barycentric subdivision and dual triangulation, and 

*^ : CP{K) C"--P{L). 
However Cp{K) and Cp{L) are both contained in Cp{BK) as we have seen. Let 

W^^ : CP{BK) ^ 0P(Af), 
denote the Whitney map. Then we have Hill, for x £ Cp{K), y G C"-p~^{L), 

1. < *^x, y >= jj^j W^'^iBx) A W^'^iBy). 
< *^y, X j^^ W^^{By) A W^'^{Bx). 

2. = ^L^K CP(K). = ^K^L C^iL). 

for the hyper-cubic complex. These are the discrete analogies of the interrelationships between d, d^, * and < •, • > in 
the continuum. 

Note K ^ L and that properties of dxjdx analogous to those for differential forms only hold if K , L are both 
regarded as elements of BK. This feature of the discretisation method is, as we shall see, crucial if we want to preserve 
topological properties of the original system. 

In the next section we derive a Whitney map in the context of a hyper-cubic complex (which matches the one 
given in |l4j). In section III, the main features of the new scheme are introduced in relation with the definition of a 
contraction operation; in section IV, we proceed with the Lie derivative, the Lie algebra and associated group action 
on the space of Whitney forms. In the last part we discuss the convergence and a numerical implementation. In the 
Appendix we outline in what sense the theory still captures topology, and display some calculations. The subsequent 
article will tackle the metric along with the covariant derivative and the curvature. A special role will be given to the 
dual complex L in relation with the metric connection. We will also propose some applications to physics. 



II. THE HYPER-CUBIC COMPLEX 



In what follows, we will consider a hyper-cubic complex, since it is the most used in the literature, and we have 
already pointed out its advantages. We consider singular homology rather than homology (see the definition of a 
simplicial complex ^3)- "^^^ describe how to handle the new setting. Consider the surface of the cube, topologically 
the two-sphere S'^. Consider the face [0123], and let us follow the steps of the construction of the Whitney map given 
in the last section. The zero-form Whitney elements are obtained by application of the boundary conditions and are 
given by: 

{Mo = (1 - - y), fii = x{l - y), H2 = xy, /is = y{l - x)}. (1) 




FIG. 2: The cube with labels used throughout the discussion. The thick loop is a holonomy around the vertex [6], discussed in 
the follow-up article. 
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The book-keeping for the co-boundary operator is as follows. Consider the simplex [0123]. Then the application of 
cIk is shown below. 

d^[0] = [10] + [30] (2) 
d^[01] = [0123] (3) 
d^[10] = [1023] = -[0123] (4) 
^-^[0123] = (5) 

Besides, the orientation is fixed by specifying a two-chain in this case [0123]. The way to check that the d operation 
is correct, is by verifying that the Stokes theorem holds. However, by using the formulae of Whitney directly given 
the cubic barycenter coordinates one finds that the requirement W^d^ = dW^ is not satisfied which then makes the 
basis list of forms an incorrect basis of Whitney elements as noted by Samik Sen . 

To remedy this problem we derive the Whitney map by induction. A Whitney map is defined as follows. Take the 
barycenter coordinates satisfying 

M.([i]) = 5,, (6) 
and generate the Whitney elements by application of 

= (7) 
<af\TyW(af)> = 5.,. (8) 

which is used as the definition of Vl^(''+i) given W'^'^\ The Whitney elements thus obtained satisfy automatically the 
Stokes theorem. 

The reader may check that the list of all Whitney elements given in the appendix A do indeed satisfy this rule (and 
may skip the rest of this section) . 

We will sketch the proof by induction, a detailed one will be given elsewhere. The first observation is that for zero 
co-chains, adopting the functional form of the example of the square above, we find that any vertex [i] has within a 
given top dimensional simplex A" of which it is a face the following expression: 

Wi[i]) = Hfc/fc (9) 

where each factor is either 

/. - /+ = (1 - (10) 

or (11) 

= /r-^- (12) 

The factors enforce the boundary conditions JHJ, and depend on which bounding simplex we are taking for the 
evaluation of the Whitney map. The claim is then that we can extract the remaining elements, i.e of higher degrees. 
An immediate observation is that application of d will suppress factors fi and replace them up to sign by a form dx^ . 
Moreover, the factor fi that is suppressed relaxes one of the boundary conditions relating to the simplex containing 
both the original one and its extension caused by application of d^ . In turn, this leads to a well-defined pull-back 

(r) 

map for integration over a given simplex which we denote by 



< a, 



. >*: Q'-iM, R)i — M (13) 



and the space of Whitney elements is generated by the image (denoted S) of 1 G M under the pull-back for each and 
every r-simplex in a^^^ G A^^\K, M), 

^MaW(X, M) 3(< A''{K), . >* 1). (14) 

This map, gives us the Whitney elements. Take a , then assume is given. Then by application of the 

pull-back map after evaluating dW{a'f'') using (0 gives for a (r + l)-simplex. This process is done until the 

list of (r -|- l)-simplices is exhausted. 
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Consider a basic example; in the appendix take /io- Then apply the exterior derivative to get a linear combination 
of terms which match VF([01]), W{[03]) in the square [0123]. Then select W^([01]). 

A side issue is that of completeness. One might be under the impression that we are calculating W for exact forms 

in 

Z'-iK, M) (15) 

Obviously, this is not true, the pullback allows us to evaluate the Whitney map on the basis set of C^{K, R) which 
is the collection of simplices A^{K). The Whitney elements obtained are the correct ones (see the appendix). 

To visualize the Whitney elements, we sketch below what a 1-co-chain looks like. As is apparent from the formulae, 
a one-chain with non-zero coefficients along the edges parallel to the y-axis have the following shape, you see that the 
form jumps over cells which are parallel to the 1-chain but this does not affect the differentiability of the form since 
d = dx A on a dy-form. and in multiple dimensions you get all derivatives apart from ^ . 



w 




X 



FIG. 3: The two dimensional plane with a square complex, coordinates {x, y) and a Whitney form with function coefficient 
sketched in the M^-£ixis. 



At this point we take the topological theory above and modify it in order to consider geometry. The key construction 
that leads us to geometry is the formulae given below for a hyper-cubic interior product and the numerous properties 
it possesses as we proceed to describe. 

III. SETTING UP THE EXTERIOR CALCULUS 

In order to extend the range of applicability of the method we need to be able to study "flows" . This involves 
introducing the idea of a vector field in the discrete setting and constructing an analogue of the Lie derivative. 
To do this we need to define an discrete analogue of the interior product. Geometrical information in the sense 
of introducing a metric in this frame corresponds to following the viclbcin formalism and is related to the Rcgge 
approach. We postpone consideration of this important construction to the follow-up article and focus on vector fields 
and the Lie derivative and algebra. The metric is the locally flat metric mentioned in section I. 

We now begin our construction by noting that the discretisation scheme is too restrictive in order to consider a 
geometry. What we need is a contraction map iy in order to consider vector fields and geometry. With d and i„ we 
will construct the Lie derivative Ly (and in turn, in the subsequent article, with some prescription for the metric and 
the Lie derivative we will be able to construct a covariant derivative V^)- 

A. Interior product 

The problem we start with is to construct a hyper-cubic analogue of the interior product (or contraction) map: 

i^w^*-) = ^^ v^uj^^^.^^Jxi^^ A ... A dx'^^ (16) 
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Our strategy is to apply the de Rham map to the dual one form of the vector field v represented as a one-chain (i.e 
a linear combination of edges) and to the degree (r — l)-form that results from lo after contraction. Thus providing 
us with a pair of chains. Let us extract it from the continuum by rewriting as 

z„w« = ^^v"" fuj,^,...^jx''^ ^...^dx^■^, (17) 

and then notice that this expression, as a pairing of two forms, can be represented on the chain space of the K complex 
as: 

where a is a one-chain, a is an r-chain and 77 an (r — 1) chain which results from contracting a. It is similar to what 
we have in the continuum with exception that a, a one-chain plays the role of a coefficient function and so should be 
treated as such. Second remark is that we will see that in fact this leads to a well-defined operation on the complex 
and is exact when the form cj and the dual to the vector field v are in fact Whitney elements. Also, we introduced a 
generalized product on the space of chains, denoted by J'^ the choice of the symbol is made to stress that it comes 
from the contraction operation. It clearly distinguishes between the chain on its right and that on its left. Let us 
postpone the formal construction and consider examples. Start with a 1 D example. Consider the edge [01] on the 
line. Then, the Whitney elements are 

{W{[Q\) = 1 - X, W{[1]) = X, W{[Ql]) ^ dx.} (19) 

and 

Z[oi][01] = [01]J^([0] + [l]). (20) 

That is, wc took the sum of the two vertices that remain when we remove the edge [01]. Also, the construction is 
guided by the embedding in the continuum: 

Winoi] [01]) = W{[01]) X W{[Q] + [1]) ^dxxl = l. (21) 

Again we said that the one chain a plays the role of a function coefficient so, in this case it is 1 and hence the result. 
If you contract [0] or [1] with the edge [01], you get zero because the edge [01] is not a face of either vertices. Consider 
the 2D example of the square. 

[3] [2] 



roi ni 

FIG. 4: Two parallel edges in a given square which arise from contraction of [0123] with one of the edges [01] or [32]. 



The Whitney elements are: 

{W^[0] = (1 - .t)(1 - y), W[l] = x{l - y), W[2] = xy, W[Z\ = (1 - x)y, W^([01]) = (1 - y)dx, (22) 
W^([12]) = xdy, VF([32]) = ydx, W{[0'i\) = (1 - x)dy, W^([0123]) = dx h dy.} (23) 

This time there are many combinations of two chains which lead to a non-zero contraction, for example the contraction 
of the square with an edge, say [01]. let us see two examples which will help us to specify iy. 



Z[oi][0123] = [01]J^([03] + [12]) 
*[oi][32] = [01]J^([3] + [2]). 



(24) 
(25) 
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And they lead to 

W^(i[oi] [0123]) = W^([01])J^1^([03] + [12]) = (1 - y)dxfdy = (1 - y)dy (26) 
W(i[oi][32]) = W^([01])J^W^([3] + [2]) = (1 - y)dx\^x = (1 - y)x. (27) 

To comment, the first calculation shows how when contracting a 2-chain, we get a linear combination of the edges 
perpendicular to the edge which specifies the vector field. This generalizes obviously to higher degree chains: we get a 
linear combination of the (r — l)-simpliccs which are orthogonal to the vector field (an edge). The second calculation 
is there because it shows that although [01] is not a face of [32], the contraction still applies because they are parallel 
within the same cell of highest degree which plays the role of the local open set. 

We formalize a little the interior product we have constructed and in the next sections, we will have to modify the 
other operators accordingly, notably the co-boundary. 

The discretised interior product is a map from the product space of chains to the product space of chains, in the 
way that iv<J is a pairing of a one-chain and an r-chain which is mapped to a priori an r — 1-chain but is a bit more 
general and the vector field is explicit: 

i : C'^^\K, M) X C'^''\K, M) — > C'^^\K, M) x C'<'-'^\K, K) (28) 
{v,cj)^vfri{a) (29) 

where rj^a) is the (r — l)-chain that results from the contraction by taking a linear combination of faces of cr^''' . There 
is also a sense in which the contraction is done in a "local" open set which is determined by the cell (or simplex) of 
highest dimension which contains both v and a. 

Also, the Whitney map applied to iyCF treats the vector field as a function coefficient i.e: 

W{v\^'n{G)) = W{v)\^W{'n{G)) = ip^{W{v)) a W{r]{a)) (30) 
where if^ deletes the constant form part of the Whitney element, 

(p^{u;dx) =u!. (31) 

The generalization of this operation when more than one contraction has been applied to a chain, like for the Jacobi 
identities below, is done in the obvious way. 

Having shown that the operation works in the cell, we return briefly to the consideration of chains, we alluded to 
them in the discussion of FIG. 3. There is a pairing of chains in i^cr which makes the map bilinear. It is clear that 
each simplex gives rise to a co-chain with support over all cells that have that simplex as a face. This is related to the 
notion of calculating in the "local" open set we alluded to. Otherwise, we are left with the consideration of constant 
vector fields, we will return to that in the discussion of the Lie derivative. We write the bi-linearity of the interior 
product as 

iycr = ^^iv,o-j X! ivi<^j- (32) 

In words, given a one simplex Vi from the one-chain v, and one term aj from the chain a, if they are both faces of 

a higher dimensional hyper-cube cr^\ then the terms contributes. It is then straightforward to establish that the 
discrete interior product squared is zero. 

Having introduced the interior product, we allow for geometry, but in doing so, the closure of the operations is on 
the product space of chains. We will establish that the other operations A, d and ★ can be modified in a way that 
makes the theory consistent topologically. Also, we will see that the Jacobi identities work exactly thus limiting the 
result of the chain operation to the product of two such spaces, while the intermediary steps involve three. So let us 
modify the discrete theory to accommodate the new operator iy. 

B. Exterior derivative 

A priori it is already under control, since if w is a Whitney form, we know how to calculate cL;, it is a linear 

combination of (r + l)-degree Whitney elements using the co-boundary d^ = on the complex. However di^cr = 
d{v\^ri{a)) is problematic, we need to make a prescription for it since it is an object in the product space of chains. 
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The outcome is that we can define it in a consistent way which also leaves the cohomology intact. To see this, go 
back to the 2D example and let 

d^(z[oi][0123]) = d^([01]J^([0] + [1])) = 2[01]J^([30] + [21]) (33) 

while mapping both sides to the continuum we get the correct result which is 

dty{{l - y)dx) = -2{1- y)dy. (34) 

The operation is well-defined in general, when applying d after i^, we only take partial derivatives with respect to the 
orthogonal directions to both w'^^-' and cr*-''-' and r > 1. So we can apply the Leibniz rule for functions in a combinatorial 
way by writing 

d^M = 2v\^ d^T^ia) (35) 

where d^ is the original co-boundary map, and thus, the homology of the theory in the product space is the same as 
that of the original theory and leads to the desirable result: 

Wid^i^a) = dW{i^'n{cr)). (36) 

To justify the formulae, we can give a proof, by taking the general Whitney form in a hyper-cube which was discussed 
above. So let 



w^(aW) = nLi/rnjr;+ j+ /\ dx'^" . (37) 

Then, starting with the right hand side of H35|l . we get 

Wi2v\''d^r,ia)) = 2(Ty(t;))J^ (^ A-dx' A)W{v{a)) (38) 

the point is that when you move the derivatives across to the left, you pick up all the terms since the vector has 
at least every fi factor that a has. There is a minor sign issue which as to be watched out for, since as stated above, 
we allow to contract with edges that are parallel to a and -^fj^ = Tl- 

Return to the co-homology. The di^j-property H35|l can be used twice to show that 

(d^)2(z„a)-?;J^(d^)2,7(a)=0 (39) 

since {d^Y = on the space of chains. We had the property that 

{d^fC'\K, R) = 0, (40) 

and we have deduced that 

(d^)^i„C"'(if, M) = 0. (41) 

The diu-property is particularly convenient since it also helps to consider exact chains as well as closed chains. Let 
us consider an (r -I- l)-chain which is exact, 

^(-+1) =rf^f^W. (42) 

then, considering a one-chain representing the vector field v and assume that we contracted an {r+ l)-chain that gave 
v\^a. Then, immediately we get 

y^K^(r+l) ^ ^ „JK (d^^M). (43) 

which is a rewriting of H42|) leading to the Homology groups 

{M, R) ^ B*'') (M, M) {M, R) . (44) 

The picture is the following: we had the link between the de Rham co-homology and the homology as a basis for 
discretising differential forms with the operator d. Now, we have extended the theory to the product space of chains in 
a way that preserves the content of the topological theory, and furthermore gives us access to the geometry. Whitney 
elements have a central role in the theory. 
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C. The wedge product and the Hodge star 

The product is much the same, but it is worth pointing out how it co-exists with iy since we have to check that it 
acts as an anti-derivation on the space of chains i.e 

/3(f)) = f3^P^ + (-l)V('') A^ iyf3^P\ (45) 

Again the criterion is given by the continuum analogue which is 

l¥(i,(aM /3(P)) = l¥(i,aW) A^ p'^P^) + (-l)''W^(a('^)) A^ W{iyf3^P^). (46) 

But we have aheady proved that because the interior product is exact on the space of Whitney forms. So by using 

i^yW = Wiy, (47) 

we rely on the continuum identity and so (|45|) . As is known, the wedge gives rise to the cohomology ring, and that 
the discrete wedge (or cup product) is non-associative, or it is only up to a factor. By acting with d on cr A^ r], 
we invoke the Liebnitz rule which is exact only after integration (application of the de Rham map) in the original 
geometric discretisation. The point is that Whitney forms were only present as dummy integration variables which 
are compatible with d, while here we want exactness before integrating thus preserving the topological properties 
associated to d and to i„. Then, to take the wedge product of two "contracted" chains is done in the natural way 
(where (g> is the tensor product) 

(CJ^ry(a)) A^ (/3J^rKa)) = (C ® /3)(77(c7) A^ v{a}) (48) 

again in a similar fashion as for d you see that we are combining the homology groups in the same way as on the 
chain space: 

H^'''> {K, R) A i? {K, K) ^ H^''+P'> (Af , M) . (49) 
as can be seen by noting that the term under brackets is the ordinary wedge operation in the chain space. 

The Hodge star has two representatives in the discrete theory, one in the K subdivision and one in the L subdivision. 
We write: 

*''{f3\''v{'j)) = f3\''{*''v{'j))- (50) 
Then, the associated properties of -k^ , -k^ and (d^)^ hold, since using H50|l . we get 

(d^)t = (★^d^*-^)(*^d^*-^) ^ ★■^(d^)2*-^ = (51) 

where we used the fact that 

■k'^*^ = IdL,*^*^ = Idx- (52) 

and in much the same way as for the sedge, we still have the Poincare duality 

(Af, M) 9i (Af , R). (53) 

What is important at this point is that the operations (A, d, k) are essentially topological and the function coefficient 
is passive or gives rise to a factor of two which does not modify the homology, a discussion of the co-homology of the 
theory in the product space is given in the appendix. Let us turn to the Lie derivative. 



IV. LIE DERIVATIVE AND CONTINUOUS SYMMETRIES 



With the operations d and iy in place, it is a matter of consistency that the discrete Lie derivative, as given by the 
Cartan formulae 

Ly = iyd + diy (54) 
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should have the following analogue: 

U = «f + d'^if . (55) 

We will find that we can go quite a long way with it and consider a discrete version of the Lie group. As we saw, we 
consider the vector field as being represented by a one-cochain which is its dual differential form. The theory is seen 
as a theory of differential forms and so this choice is imposed to us. If we were to consider the vector fields in the 
dual complex as -k^v in the L, we would loose the common support that both the vector field an the form must have 
in order to carry through the construction of the exterior calculus. We start with examples, since we already have 
the operators. 

Example 1: (see the Appendix), to investigate the construction, we start with the vector field (dual to) H169|l and 
the form ^if^ . These are the continuum analogies of H175|l . H176|l . The vector and the form will be considered on 
three faces, [0176], [1247] and [0653] for this purpose. To each of those corresponds a volume element two-form, as 
given in section II (i.e a constant form such as dx A dz). Note that due to the orthogonality of the form with the 
vector chosen, we get iyV^ = 0. Next, evaluation of dv^ in the discrete setting H178() is obtained by means of applying 
the rules for the hyper-cubic complex given in section II. After application of the Whitney map the result (e.g H181|) ') 
is obtained by application of the rule: 

This means that the transcription into an ordinary differential form is obtained by suppressing the form part (like 
a da; or a dx A dy) on the left hand side of the J ^ symbol. So by substituting the list of Whitney elements given in 
Appendix B, and the rule, we obtain H181|) which leads to the identification of the various terms in (|173|l . We have 
found the following: 

W{Uvb)=LwvW{vb). (56) 

We now move to a more delicate case, for which we will need to clarify the notion of the tangent space TuM of 
vector fields, where U may be taken as an open set taken to be the cell. We wish to point out that we do not think 
of the tangent space TM as a fiber bundle. The topology of the complex is too restrictive for that, open sets only 
overlap on a common boundary. However, the Whitney elements, because their support overlaps over boundaries 
guarantee as we will see, that we can consider properties of vector fields. 

Example 2: Take the vector u to be parallel and in fact identical to . In this case, one has to be careful since 
if we use the proposed formulae for the interior product, we get (|188|) . which is clearly wrong since there should be 
some terms with coefficient vz-^^vz^ as shown in the appendix. This is the point at which we introduce the rule above 
from which we get (|189|l . Considering (|189|l . keeping only the terms in [0617], 

iuVb - (4. - vz,vzM\''m + [6]) + (4, - vz,vz,)[n]f{[l] + [7]), (57) 

then, 

diuVb = 2(4^ - vz^vz^mW^'im + [76]) + 2(4^ - vz,vz,)[n]fm] + [67]) (58) 

= 2(-(42 - vz,vz2){l -x) + (wli - vziVz2)x)dx (59) 

= 2{vzi - wz2)(wz2(l - x) +vzix)dx (60) 

as required to match the continuum expression. For the i„(i-term, the problem of support which lead us to modify 
the coefficients does not arise and we get 

ludvb - (-4^ [06] - vz.vz, [17] + 4^ [17] + vz.vz, [06])J^([01] + [67]) (61) 
= (-42(1 - x) - vz^vz^x + v%^x + vz^vz^i^ - x))dx (62) 
^ [vzi - vz2){vz2{'^ - x) + vzix)dx (63) 

which gives the correct matching on [0617], not forgetting the multiplicative factor of two specified by the rule for 
"div" . This gives us the correct matching. 

The combination of Example 1 and Example 2 cover all cases of a pair of a form and a vector field that can be 
expressed in the present framework (i.e parallel or orthogonal) . They constitute the two parts of the proof that the 
discretised Lie is the correct one. What we are describing here is a flow from the face [0176] to the face [0653] and 
check its vertical (z-component) at the edge [06] which joins the two faces. The tangent vector to the integral curve 
is thus rotated by a right-angle as it propagates along its flow line. The component lives in [0176], the component 
Vy lives in [0653] while the component Vz lives in both. 
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Evaluating components of the discrete Lie derivative as we have done above will be particularly useful when we 
turn to discuss the covariant derivative in the subsequent article. But before we address this issue, we ought to check 
how the Lie algebra property, and notably how the Jacobi identities translate to the lattice. This property gives us 
a lifting of the Lie algebra of vector fields (which we have yet to define) to the algebra of derivations acting on the 
space of forms. In the continuum we write: 

[Lv,Lu]=Ly^y (64) 

Although we are in a position to evaluate the left-hand side of (j64|l . we need to make a prescription for the complex 
based bracket [v, u], and in fact we know that in the infinite dimensional representation of vector fields which we 
adopt, the continuum commutator is clear: 

[u, v] = {u'd.v^ - v'd^u^)dj. (65) 

Or conversely, we need to show that starting with [L^,, L„] , the vector field which is the argument of the Lie 
derivative on the right-hand side is indeed the discrete analogue of [v, u] given in (|65|l . 

Given two one-chains a and 77 on a 2D complex, we introduce their bracket 

[. , .]^ : C\K) X C\K) C\K) x C\K) (66) 
(M, [77])k-> (67) 

where in components, assuming summation over the indices, we are given a = cr^ [zj] and 77 = r]ij[ij], then, 

[v, u]^ = {cTijilki - nij<^ki)[ijkl]\'^[ij]. (68) 

Example: Again, much is achieved by an example. Take the face [0671], and using the previous example, we take: 

ii = vx,m + vx.W] (69) 

v = vz,m+vz^[l7] (70) 

The continuum analogies are 

Wu — vxi (1 — z)dx + vx2zdx (71) 

Wv — vz^i^ — x)dz vzixdz (72) 

Application of 165|) . leads to 

\Wu, Wv\ = (wxi(l - z) + vx2z){vzi - vz2)dz + {vz^i^ - x) + vzix){vx2 - vxi)dx (73) 

It is straightforward to reproduce this formulae on the lattice by using the prescription. 

The construction can be clarified further by considering the example of so(3), the rotations in 3D. We will setup the 
Lie algebra and show that it is well-defined on the complex. Second, we will explain how, although the Jacobi identities 
are found to hold, the lifting to the group is only approximate which is because integral curves are approximated by 
Whitney elements. Consider the so(3) generators and the form lu respectively given by 

X^yd,- zdy, (74) 

Y = zd,-xd,, (75) 

Z = yd^-xdy, (76) 

UJ = Wi[65]) = (1 - x)zdy. (77) 

A direct calculation in the continuum leads to 

L[x,Y]Wi[65]) - [Lx, LY]Wi[65]) = -z(l ~ x)dx - yzdy (78) 

which are the continuum Jacobi identities for the Lie. We now proceed to analyse the discrete version of the Jacobi 
identities by constructing the complex based version of this calculation. We will use the rules we have introduced, 
notably the c?i„-rule, the rule for immersion and the rule for products of coefficients given below. The central feature 
is that of closure, although we use a product space, we can recover the exact expressions by simply identifying: 

1A4>^4>, (79) 
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after applying the latter rule. In the 3D coordinates available, the generators are exactly represented by: 

Y = ([17] + [24]) - ([67] + [54]), X = ([35] + [24]) - ([65] + [74]). (80) 

because a look at the table of Whitney elements in the appendix, will help to see that these are the dual one forms 
of the vector fields after mapping them to the chains to the continuum. The bracket we proposed was guided by 
requiring 1) The closure of the algebra under the bracket, 2) The functional form of [X, F], 3) The Jacobi identities. 

[X, [y, z]] + [y, [z, X]] + [z, [x, y]] = o (8i) 

and at the level of the Lie derivative it reads: 

L[X,Y] = LxLy — LyLx, (82) 

we will derive this relation below for the complex based theory. This means that we seek a complex based bracket 
operation that maps to the continuum one by a Lie algebra morphism. With the basis of forms we have, we can 
already set what the answer is: 

[X, Yf = AiiWiX), W{Y)]) (83) 

and this is satisfactory because it is an exact relation and not an approximation as: 

W[X, Yf = Z (84) 

the continuum Z, which means that W is a Lie algebra morphism. So this would be enough to establish 2), and 
certainly proves the existence of the discretised bracket. 

What we have is that for a Lie algebra g in the infinite dimensional representation of vector fields. The map W is 
a Lie algebra morphism if and only if, the collection of one form dual to the generators of g can be written as a linear 
combination of Whitney elements. 

The bracket in H83(l is an implicit definition making use of the continuum one. It is satisfactory as we were able to 
give an exact analogue of the vector fields X, Y, Z. We proceed to the calculation of the Jacobi identities for the Lie, 
and first of all we need the list of edges that are mutually parallel, so that we will know when to apply the rule called 
"Tangent space of vector fields" : 

Li = {[17], [06], [35], [24]}, = {[67], [54], [32], [01]}, L3 = {[65], [74], [12], [03]}. (85) 

The detailed calculation is done in Appendix C, since all the rules we introduced play a part, we displayed in curly 
brackets the associated continuum calculation for the purpose of illustration. The picture should then be clear, the 
algebra of vector fields is local and is exactly captured. The Jacobi identities for the algebra 

ad{[x, y]) = [ad{x), ad{y)] (86) 

which are satisfied here, have a "local" analogue in the Lie group, consider the vector fields X^, Xy, satisfying 

[Xx, Xy] = X^x y] (87) 

and the flow of X^ is given by 

$a;(i, a) = aexp(tx). (88) 

But the latter flow is only approximate since it cannot be captured in the space of Whitney elements (for example 
a vector field rotating in the x, y plane. The limitation is easy to understand, owing to the exactness of the vector 
fields in the theory, the algebra is exact. However the space of co-chains as being given by the Whitney elements does 
not lead to the correct integral curves and the associated relations, we only obtain an approximation, which is for the 
flow of a vector field and its discrete counterpart <tx is controlled by correction terms of order — W{'tx). 

In conclusion to the part on the Lie algebra and lie groups, we gather the results. A Lie group action on the space 
of chains and associated Lie algebra of vector fields are such that: 

1) The Lie algebra is exactly captured if the generators are elements of the space of Whitney elements. 

2) The associated Lie group G is exactly captured if the integral curves of the vector fields leave the space of Whitney 
elements invariant. Otherwise it is approximated by the Whitney forms. 

We now address the central issue of implementing the technique we have just presented and expose some numerical 
results. 
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V. CONVERGENCE OF THE MODEL AND NUMERICAL IMPLEMENTATION. 



We have seen that the results are exact on the space we have introduced, however the Jacobi identity for the Lie 
is only satisfied in the local picture, flows in the group are only approximate in general. Overall, this means that we 
have good control over the convergence. 

An issue related to convergence is that of refining the triangulation (which is here made of squares). A discussion 
of this in the context of simplicial categories J^lj is worth considering. In that picture as here, one works with chain 
complexes C. and their dual co-chain complexes and it is argued that they are to be preferred to homology. Due to 
the Whitehead theorem they are invariant for connected spaces. That is, we take subdivisions of the original complex 
and one can relate the associated chain spaces. In turn, given the chain space and the dual cochain space, one can 
define the various mappings we have introduced above. 
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In the present perspective, we insist on the fact that refining the triangulation is the only mean by which one can 
obtain a better approximation to the various forms which we always take to be Whitney elements. 
Next, we discuss the example of a one-form field on the torus, satisfying the equation 

(L„ - riA)uj = (89) 

which is relevant to fluid dynamics [2L |. and is a flrst step toward setting up the Navier-Stokes equation in this 
formalism. We start by tackling the problem analytically and then numerically. The square "triangulation" of the 
torus is shown above: 

Now, the vector field is constant, taken along the a;-axis in the plane of the page v — Vxdx- We seek a one- form 
solution as a linear combination of Whitney elements. The Lie derivative is then: 

L^uj = {y+ - y^)vxdy (90) 

The notation we use for components, is such that it can be applied to every 2D cell. The sign and i or j indices 
denote a specific edge. For example in [56, 10, 9], is the component on the edge [9, 10] and is the component 
on the edge [59]. Similarly when a component outside the cell is needed, we take for example y^j to be the [10, 14] 
component. Next is the Laplacian term: 

dd^uj = (-(x+ - a;+ _j + y_ ^ - ?/_) + (a;+,, - x+ + y+j - y+))[x+] (91) 

+ - + j/_ J - y^) ~ {x^ - + - 2/_ _j))[?/-] (92) 

+ ((a;+,» - a;+ y+j - y+) - {x-^^ -x^+y+~ y+ _j))[y+] (93) 

+ {{x-.i -x-+y+- y+.-o) - {x- - x--i +y- - y--j))[x-] (94) 



and 



dUuj = ((y+j - - y_j + x+) - (y+ - x+ - y- + X-))[x+] (95) 

+ {{y+ ~ x+ -y- + X-) - {y- - x+-t - y_._j + x_ _,))[y-] (96) 

+ {{y+,t - x+,t - y+ + x-^i) - {y+ -x+- y- + x-))[y+] (97) 

+ {{y+ -x+-y-+ x+) - {y+-j -X- - y--j + x--i))[x-] (98) 
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The first solution is the constant zero-mode solution. Then, is zero because the y-edges have equal coefficients 
across each ribbon. Which means that: 

y--i = y-=y+ = y+,i = Y (99) 

y-j = = y+ (100) 

y-,-3 = y+..-, = y- (101) 

In this case, equation (|89(l reduces to: 

Aw = 0. (102) 

To get zero we need either each component be zero or each pair of coefficients (e.g the two components) to 
cancel out. Set all x to be zero and we are left only with: 

l\uj = (y+ -2Y + Y^)[y-] + {Y+ - 2Y + Y-)[y+] (103) 

Our equation is finally 

y+ - 2y + y_ = (i04) 

This gives us a constant form ui — C dy as an exact solution. Let us think about H104[l . After all this is an approximation 
to the second derivative: 



{d^fL0y = Y+-2Y + Y^ (105) 
by using a finite difference rule for the derivative: 

D^u ^huj{x + h)-oj(x)) (106) 
h 

which is correct to order h. Since [y~] has the same coefficient as [y+\ makes the associated form one with a constant 
coefficient and so having Whitney elements presents no benefit from the usual finite difference approach in terms of 
approximation in this case. 

In order to move to a more general example, we consider again a solution in y-components, but this time we let y 
increase as we move from one y-ring to the next one, that is in every cell we let y+ — y- = 0. Then, 

LyUj = {y+ - y^)vxdy = Ovxdy (107) 

In order to cancel the Lie term, we require as before that all x-coefficients be zero and so 

(108) 
(109) 
(110) 
(111) 

Equate H92|l and H93|l . We find that they have equal coefficients after noting that using H108|) we can rewrite H92|l as 

y-j -y- -y-+ y--j = y+j -y+-y-+ y-.-j (112) 

similarly for (jHSJ- Next, equate (EHl and (P7|) . 

y+ - y- - {y- - y--i) = -y+-y+ + y- (ii3) 
y+-y- = \{y+,^-v-,-^) (ii4) 

We are now in a position to write down a new solution. Let y be constant along each ring parallel to the y-axis. Also, 
we have seen that 

y+-y-^e (115) 



y+ 


- y- 


= y+.j - 


y-,3 


y+ 


- y- 


= y+.-j 


- y-,-j 


y+ 


- y- 


= y+.j - 


y-,j 


y+ 


- y- 


= y+.-3 


- y-.-j 
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If wc enforce this for all squares we have a local (in the sense of a single simplex) matching of Whitney elements which 
is fine. However, if we are to look at boundary conditions, particularly for closed surfaces like the torus, we will need 
some more prescription to guarantee closure of the field. To this end we consider 

(l) = ae'%dy (116) 

where a and b are constants to be determined. Then, in the interior of a given square, we have 

(j) = a{l + i^)dy + 0{x'') (117) 
and letting x = ix, leads to the usual expression, after we enforce the following better ansatz: 

<^ = </)(Zi)(l-|) + </.(?2)(|) (118) 

(t>{li)= j ae'^dy = ae'^. (119) 
J\h\ 

Moreover \ih = h the edge length of a small square, then we can match the function with the Whitney map applied 
to a chain field (p. We the have the matching 

<j> = W{^) + 0{x'^) (120) 

To move with the analogy one step further, note that one retains a bijection between the function e'^ and ix. i.e, the 
function cf) oscillates over a large number of squares: 

||y+-y_|| «27r. (121) 

To guarantee closure, let be the number of edges along a circle of the torus (parallel to the x-axis), then require 
closure (recall y+ — y- = 6 a constant): 

{y+-y-)N, = 2kn. {fc = 0, 1, 2, . . .} (122) 

Then, fc = leads to the zero mode we already got, the rest give the k-th mode solution. Let us do the contimmm 
calculation: The change of variable x ix leads to dx idx, ^ ~*^) * —iv^- Apply this to the solution: 

Ly(j) = vj--^e^^ dy = Va:^e'^ dy (123) 
oix b 

■kd-kd(t>= ^e'T^ldy) (124) 
b'^ 

This solves the equation provided ^ = rj. In order to move closer our analogy between the continuum and the 
complex we need the following, evaluate: 

dln(j)dy = ^dx Ady (125) 

-kd\n(t)dy = ^ (126) 

iydln (j) dy = dy (127) 



Then, 



d*dln(j)dy=^-^—^-^dx (128) 



^d^d\ncj)dy= ^-i—^{-dy) = ^(-dy) + ( ^ ) dy. (129) 



Note that according to the discrete theory we get: 

0/ gixfc+i _ gixfc 



= i{xk+i - Xk) + 0{{xk+\ - Xkf) (130) 
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Using the modulus, 



Therefore, 



l2 2 



1 



In the discrete approach this is written as 



A\ln(l>) = > / '-dy[jij2] = > ^[jij2] 



Then apply the discrete version of the operator 



(131) 

(132) 
(133) 

(134) 

(135) 



and solve the matrix equation. Next, with a varying vector field, the solution is modified, that is the rate of change 
of the solution starts varying. This does not modify the calculation above with the minor exception that now 



< sup 



' f{Xk) ' 



, 1 



0(e2) ^ 0(e2), 



(136) 



which leads to the same equation as before with a variable v^. 

Regarding the inclusion of modes, we have studied the function In (j) where in principle we could include all the 
modes (along the x-circles): 



Let 



Then we will get 



(A-L,)ln(x) = ^(-x"-^^xX') 



1 

X L 



Substituting functions, note the term for f/.. 



fk = a + hx + cx^, 



leads to the following equations for the constant term and the x-term: 

an^ + (2& — uoa)K — 2c — v^b = 
bn^ + (4c — vob — Vxa)K — Vxb — 2vq = 



(137) 

(138) 

(139) 
(140) 



(141) 
(142) 



(143) 
(144) 



These equations are to be solved for a, b and c in principle. We can then substitute back the solution into our trial 
solution: 



X = {a + bx + cx )e' 



(145) 
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thus 

X' = {nia + bx + cx"^) + b + 2cx)e'^'' (146) 



The X term in this expression is then 
The approximation is then given by: 



(6 + 2c + aK) (147) 



X = (< / >o + < / >i x)e^- (148) 

=</>o+(</>oK+</>i)a; + 0(a;2) (149) 

So 

X' = (< / >o < / >i)(a:„+i - Xn) + 0(e2) (150) 

which is to be compared with H130() . The observation being that the theory captures the hnear term. We leave here 
such considerations of matching the continuum solution. 

To gather our findings we note that the key features are the use of the logarithm and the selection of the mode 
by specifying the relevant boundary condition. The k-th mode is given by specifying the constant C which is a 
normalisation and setting: 

<j>+ = C (151) 
(t)-^C + 27Tk (152) 

Numerical tests have been done on the torus we have displayed above (with a few extra squares only), and although 
we have already discussed the convergence, we would like to explain the numerical procedure. It consists of using a 
singular value decomposition (SVD) algorithm [2^ . 



One form field solution with values in the range [400, 1400] over circles parallel to the x-axis. The solution is 
constant over y-circles, e.g [1], [7], [13], [19] are the vertices of one such circle. 

Two key advantages to this method are first its well-known stability and second that it solves a matrix equation by 
minimizing the square norm below, so in some cases were the present theory cannot be solved exactly as a matrix 
equation, the algorithm gives us a good approximation. The SVD is used to solve 

Ax = b (153) 

and leads to the following decomposition (the Wj are the singular values): 

A-^ ^ V.[diag —]U^ . (154) 

Wj 
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The solution is then obtained by minimizing 

X^ = \\Ax-bf. (155) 

We also note that the formal errors due to the fit, let V(i) . . . V(m) be the columns of V then the boundaries of the 
confidence ellipsoids are given by: 

Ax' = wliV^.ySaf + ... + wl,{V^M)-Saf- (156) 

In the case of the graph above, we took the following vector field: 

=4, Vx e [0, 3[, Vy; (157) 
v, = -2,yx€[3,6[,yy; (158) 
Vy = 0,yix,y). (159) 

Note that the change in the vector field results in a change of concavity in the graph as the boundary conditions are 
enforced on the solution. Here C = 400 and we identify = 400 with 0~ = 1400. In this example the matrix A is 
96 X 96. For a constant vector field, we just get a straight line in each section. The inclusion of a non-zero mode is 
achieved by modifying the matrix in a way that take into account the qiiotienting which lead to the principal value 
of i4>. The modes are then treated one at a time, by specifying the appropriate matrix. 



VI. CONCLUSION 



Wc have proposed a new discretisation scheme for differential geometry which is deeply related to the corresponding 
continuum theory. The first step which is central, was to accommodate an interior product. This required the 
introduction of a product space. Then, preserving the de Rham complex motivated us to consider a ring of operators 
composed of i^,, d and and the introduction of a special rule for diy, for which a prescription was required. The 
Lie derivative, as an element of the ring, required us to modify the coefficient of the chain so as to accommodate the 
case in which the vector field varies over the cell as expressed by the sum of Whitney elements. The bracket was 
also introduced for vector fields in the infinite dimensional representation, for which the geometrical interpretation is 
crucial. Also, the Jacobi identities for the Lie is exact, but one bears in mind that in general the flows are approximated 
by Whitney forms. We also have some numerical implementation for the torus and a hold on convergence formally. 
We are now in a position to introduce metrics in this scheme which we proceed to do in the article 11. 



VII. APPENDIX 



A. Hyper-cubic Whitney elements 



For the purpose of this article we need the list of Whitney elements given in the table. 



Zero-cochains 


One-cochains 


Two-cochains 


^o = {l-x)il-y){l-z) 
fxi =x{l-y){l-x) 

fi2 = xy{l - z) 
M3 = y{l-x){l - z) 
fj,4 = xyz 
fir, = (1 - x)yz 
Ha = {1 - x){l - y)z 
tn = x{l - y)z 


VK([01]) = dx(l-j/)(l-2) 
W{m) = dy{l-x){l~z) 
iy([06]) = dz{l-x){l-y) 
W{[n]) = x{l-y)dz 
W{[12]) = x{l - z)dy 

W{[2A\) = xydz 
W{[2,2]) = 2/(1 - z)dx 
1V([35]) = 2/(1 - x)dz 
Ty([65]) = z{l - x)dy 
W([67]) = (1 - y)zdx 

Vy([74]) = xzdy 
W{[hA]) = (1 - z)zdx 


W([0123]) = (1 - z)dxAdy 

W([1247]) = xdyAdz 
W([0356]) = {l-x)dxAdz 
Vy([0167]) = il-y)dxAdz 
W{[6745]) =zdxAdy 
W^([3245]) =ydxAdz 
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B. Relation to the de Rham complex 



Let us discuss the relation of the setup to the de Rham complex, which is extracted from differential forms under 
the exterior derivative d. Having shown that on the extended space (after application of iy), (d^)^ = is a priori not 
sufficient. But note that we have apparently lost one important feature, which is the de Rham theorem which 
states that the map A, defined as 

A : HriM) X H''{M) — > K (160) 
is bilinear and non-degenerate. This map could be expressed as 

(a, W{r,)) = f Wiv). (161) 



However, with the construction introduced above, we seem to have lost the non-degeneracy property of A. To see 
this, let C G CiK, R), and let tj^''+^^ G C+^iK, R). Then using the rule for immersion, we get 

/ ip°{W{v)) AW{ii,r]) (162) 

J a'' 

to be compared with 

WiC). (163) 



Clearly, the space of forms spanned by products of chains is larger, but we have shown, in the proof of the Hyper-cubic 
rule that the action of d^ on the larger space has some functorial properties which link it with d, we found that, 

Wd^iy = dWiy (164) 

and so we may put the various steps of that calculation under the integral symbol. It means that we have a larger 
space of forms and that the co-homology is "intact" . This is another virtue that we extract from the din-rule. 



n^'iM, R) ^ (Af, R) 

"a |a (165) 

CP {K, R) C {K, R) ® CP {K, R) 

A comment on the diagram: we wrote A where we could have written W, but we stress that we are using the inner 
product of a space of chains with the space of differential forms. The vertical arrow on the right relates two larger 
spaces than that on the left. The diagram being commutative entails the compatibility of the present scheme with 
the original geometric discretisation. One could say that the introduction of geometry does not modify topology. 

In terms of the homology we got, 

H''{K X K,R, d^) ^ H''{K,R, d^), (166) 

while, 

B'^iK, R, d^) <B''{Kx K, R, d^), (167) 
Z''{K, M, d^) < Z''{K X K, R, d^). (168) 

Hence, the extended space of chains consisting of elements of the form specified by the rule for immersion form 
a larger subspace of the space of forms under the modified Whitney map than does the ordinary chains under the 
Whitney map. The extended space can be embedded into the product space which has the same topology. 



C. Examples 

Example 1: Let 



vxi (1 — z)dx + vx2zdx + (1 — z)dy + iiy-^ziij/ (169) 
^22(1 — x)dz -\- wz2(l ~ y)dz -(- vzixdz + vz^ydz (170) 



23 



Note that 



iyV^ = 0. 



Next, 

dv^ = —vz^dx Adz — vz^dy Adz + vz^dx Adz + vz^dy A dz, 

Lyv'' = iydv'' = -vz2 {vxi {1- z) + Vx2z)dz - vz^ {vy-, [l- z) + VY2z)dz 
+ VZi{vXi{l - z) + VX2Z)dz + vz^ivYiil - z) + VY2Z)dz 

Lattice calculation: 

V = vx, [01] + VX2 [67] + VY, [03] + vy^ [65] 
Vb = vz2[06] + vz^[17]+vzs[i5] 

Again, see the discussion of the interior product, 

iyVb = 

Then, 

dvb = {vz, - vzj[6017] + {vz2 - vz3)[0653] + vz, [1724] + vz, [3542] 

and 

t,{vz, - vz^neOU] = {vz,-vz2){vxAm\''m + [01]J^[17]) + z;x.([67]J^[06] + [67]J^[17])) 
iv{vz2 - vz,)[0653] = {vz2 - t^Z3)Ki([03]J^[06] + [03]J^[35]) + t;y,([65]J^[06] + [65]J^[35])) 

the continuum analogue of which is obtained by identifying the sum of Whitney elements with a constant form 

Wii,{vzi -wz2)[6017] = (wzi - vz2)ivxi{{i - z)dx\dz) + vx2{zdx\dz)) 
Wii,{vz2 - wz3)[0653] = -{vz.2 - uz3)uyi(((l - z)dy\dz) - wy^ (zdyjdz)) 

Which matches the continuum expression after use of the continuum identification, simply by deleting the J sigi 
the constant differential which is part of the coefficient form. 
Example 2: Next, try with u which is parallel to Vb- Then, 

iuVb = {vz2{'^ -x) + vzixf + (t;z2(l - y) + vz^yf 

and, 

diuVb = 2{vzi - Vz2){vz2{^ - x) + vzix)dx + 2{vz3 - vz2){vz2{^ - y) + vzsy)dy 
Next dvb was given above and, 

iudvb = {vz2 (1 - a;) + vz^x){vz2 - vzi)dx + {vz2 (1 - y) + Vz3y){vz2 - vz3)dy 

and finally, 

LuVb = {vzi - yz2)(.uz2{^ -x) + vzix)dx + {vz-^ - wz2)(^'Z2(l -y)+ vz2y)dy 

Lattice calculation: 

u = vz2 [06] + vz^ [17] + vzs [35] = Vb 

Then, 

iuVb = 42[06]J^([0] + [6]) +4J17]J^([1] + [7]) +43[35]J^([3] + [5]) 



171) 

172) 

173) 
174) 



175) 
176) 



177) 

178) 

179) 
180) 

181) 
182) 

and 

183) 
184) 
185) 
186) 
187) 
188) 



in order to reproduce the cross terms we set: 
and 

d^hv, = (4J06] + i;zi«zJ17])J^([10] + [76]) + (4, [17] + «z,«zj06])J^([01] + [67]) 
Next, we aheady calculated dvb, so that after retaining the components in [0653] and [0176], 

2iiidvb = -vzAvz.m + ^z,m)f{m + m) + vzA^z^m - vz,m)f{m + m) 
+vzAvz,m+vz,m)\''m] + m)+vz,{{vzsm - vzA^m^'m + [65]) 

Gathering the terms, we get 

iiidvb = {-v%,m-vz,vz,m+vi^[n]+vz,vz,m)\''{M + m) 

+ (4J06] - vz,vz,m + vl,m - ^;z3^;zJ06])J^([03] + [65]) 
Example 3: The complex based calculation is as follows: 

»^[65]=-([65] + [74])J([6] + [5]) 

iy[65] =0 

d^ix[&f>\ = -2([65] + [74])J([06] + [76] + [35] + [45]) 

{ = -2z(l-x)dz + z'^dx} 
d^[65] = [6574] + [6530] 

{ = (1 — x)dz A dy ~ zdx A dy} 
i^d^[65] = -([65] + [74])J([67] + [54]) + ([65] + [74])J([06] + [35]) + ([35] + [24])J([65] + [03]) 

{ = y{l — x)dy + z{l — x)dz — z^dx} 
i^rf^[65] = (-[67] - [54])J([65] + [74]) - ([17] + [24])J([03] + [65]) 

{ = —x{l — x)dy — z^dy} 

So, 

Lx[65] = -([65] + [74])J([06] + [35]) + ([35] + [24])J([65] + [03]) 

{= y{l — x)dy — z(l — x)dz} 
i:y[65] =iyd^[65] 

Next, 

iyLxm = ([65] + [74]) ® ([17] + [24])J([0] + [6] + [3] + [5]) 

+ ([65] + [74]) ® (-[67] - [54])J(-[7] - [6] - [4] - [5]) 
{ = 2^(1 — x)x — xz{l — x)} 
d^iYLx[65] = 2([65] + [74]) O ([17] + [24])J ([10] + [76] + [23] + [45]) 

+ 2([65] + [74]) ® (-[67] - [54])J(-[17] - [06] - [24] - [35]) 
{ = 2^(1 — x)dx + 3z^dz — 2zxdx + 2(1 — x)xdz — z{l — x)x + xzdx — 3z^dz — x{l — x)dz} 
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then 



d^ixd^[65] = ([65] + [74])J ([6701] + [5423]) - ([65] + [74])J ([0617] + [3542]) (216) 

- ([35] + [24])J ([6574] + [0312]) (217) 
{ = —zdx A dz — 2zdz A dx — ydx A dy} (218) 

iyd'^ixd^lGb] = 2([65] + [74]) ® ([17] + [24])J (-[01] - [67] + [54] + [32]) (219) 

+ 2([65] + [74]) (g) (-[67] - [54])J(-[35] - [42] - [06] - [17]) (220) 

+ (-[35] - [24]) ® (-[67] - [54])J ([65] + [74] + [03] + [12]) (221) 

{ — —z'^dz — xzdx + 2z^dz + 2xzdx — yzdy} (222) 

^^^^^"^[65] = ([65] + [74]) (-[67] - [45])J(-[6] - [7] - [5] - [4]) (223) 

- ([65] + [74]) ([17] + [24])J([0] + [6] + [3] + [5]) (224) 
{ ^ ~xz{l ~ x) - z^} (225) 

d'^iYixd^lGb] = ([65] + [74]) (-[67] - [45])J(-[06] - [17] - [35] - [24]) (226) 

- ([65] + [74]) ([17] + [24])J ([10] + [76] + [23] + [45]) (227) 
{ ~ —z{l — x)dx + xzdx — iz'^dz — x{\ — x)dz} (228) 



LyLx[65] = — z(l — x)dx + zxdx — iz^dz — x{l — x)dz (229) 

+ 2z(l - x)dx + Sz^dz - 2zxdx + 2(1 - x)xdz (230) 

— z^dz — xzdx + 2z^dz + 2xzdx — yzdy (231) 

= z(l — x)dx + x{l — x)dz + z'^dz — yzdy (232) 

Exact matching is also found in the calculation of LyLx[65]. 
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